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Abstract 

Since nearly 4 decades, various theoretical behaviours have been found for the ther- 
mopower in the variable range hopping regime. In 1969, Cutler and Mott have pre- 
dicted a linear variation with temperature T of the thermopower: S = const. T. In the 
seventies, it has been found by Zvyagin, Overhof and Mott that S = const. T(l/2). 
In 1986, Triberis and Friedman have found S = const. T(-l/4) . But there is up to 
now no theoretical formulation of the thermopower when this one is T-independent. 
By choosing a specific distribution for the density of states, we show in this paper 
that all behaviours above can be unified in a unique thermopower formula. We find 
in addition with this formula, a T-independent expression given by: S=(L/xi)(k/e), 
in which xi is the wave function decay length and L is a characteristic length, de- 
pending on the form of the density of states. 
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1 Introduction 



Mott's paperl [T] of 1969 has played a key role in tlie construction of tlie im- 
portant theoretical edifice [2]-[n] on which is rested our understanding of the 
hopping conduction mechanism. The electric properties of various disordered 
systems and amorphous semiconductors have been explained |lOj during more 
than 3 decades by using Mott's variable range hopping (VRH) theory. Such a 
theory is still extensively used today in varied fields, to investigate the mech- 
anism of charge transport for example as well in correlated electron systems 
pT] as in biological systems as DNA |12j . 
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However, the study in the VRH regime of some thermoelectric properties 
such as the thermopower have been somewhat neglected. It seems indeed, 
that this one is ill adapted in its present form to support the experimental 
measurements of the Seebeck coefficient carried out on different VRH systems. 
A simple illustration of this is its persistent disobedience during more than 3 
decades, to the predicted theoretical laws (l]-[8], observed in most amorphous 
semiconductors [6], [8], [13]- [17]. 

Because of new environmental problems, the field of the thermoelectric power 
generation has attracted recently a major interest in materials with high ther- 
moelectric properties [TS]. In this context, the explanation of the experimental 
behaviour of the VRH thermopower remains still problematic since there is 
not yet a clear theory that explains some TEP behaviours when conduction is 
by variable range hopping, as the divergent one observed at low temperature 
in ref [19] or the temperature independent one, reported by the authors of ref 

m- 

Such a theoretical insufficiency of the TEP is attributed in this work to the 
fragility of the original hypothesis concerning the nature of the density of states 
(DOS), with which the evaluation of the Seebeck coefficient has been made. 
By reconsidering the hypothesis of the "slow linear variation of the DOS" 
near Fermi level on which most of the VRH theories of the TEP are rested 
[l]-[8], we investigate by using the percolation theory both conductivity and 
thermoelectric power when the DOS takes the asymmetric generalized form: 



N{E) = N{EF)+g{E) (1) 

where N(Ei?) is a finite density of states at Fermi level Ep and g(E) represents 
the asymmetrical part of N(E): 



g{-E) = -g{E) (2) 



This choice of g(E), is motivated by the following reason: by using the concept 
of the random network of conductances [2],[S]5[H]5 we can expect to find Mott 
T^^/^ conductivity again since the number of conductances that belong to the 
critical path [3], [6], generated by N(E) of eq.([T] -121) would be quasi similar to the 
one generated by a constant DOS |2]. If so, we can expect to obtain different 
thermopower classes that correspond to the Mott T~^/^ conductivity, rather 
than the unique an insufficient formula obtained by the classic VRH theories. 
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2 Percolation method for the VRH problem 



It has been shown by using the percolation theory [2]-[Z|, that the conduc- 
tion problem between localised states is equivalent to the conduction problem 
through a random network of conductances. Each conductance aij defined by: 



, 2rij \Ei\ + \Ej\ + \Ei — Ej\ 
an = expj — - — ■ — \ (3) 



links two sites located at energies Ej and Ej, separated in space by the distance 
Tij. In eq.(I2]), k is the Boltzmann constant, T is the temperature and ^ is the 
decay length of the wave function. 

To solve a conduction-percolation problem, two critical conditions must be 
satisfied: the first concerns the number m(Ei) of incoming conductances (or 
bonds) to a site located at energy Ej and positioned at the center of a sphere 
of radius r,.,- : 



Ayr r 

m{E,) = —j rlN{E,)dE, (4) 

When averaged by a weighting probability function, m(Ej) must reach the 
critical concentration c of conductances per site, solicited for conduction within 
an active energy layer A, given by (percolation criterion): 

A 

j m\Ei)N{Ei)dEi 



I m{E,)N{E,)dE, 

-A 

The second critical condition concerns the nature of conductances that are 
solicited for conduction. To be an efficient conductor, each conductance aij 
must be larger than a critical conductance given by: 



ac = exp{- — } (6) 

Physically, when these two critical conditions are satisfied, it appears a con- 
tinuous (critical) path of conductances, joining one side of the VRH system 
to the other. The problem of the random network is said solved only when Cc 
is identified, and this occurs when eq. ([5]) is solved with respect to A. 
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3 Conductivity of the asymmetric DOS 



We show now that our prediction to find again Mott T"^/*^ conductivity for 
the DOS given by eq.([T],[2]) is true. We take for this a generalized asymmetric 
DOS form given by: 



N{E) = N{Ef) + Sg.E'^ = N{Ef) [1 + i^g.E^] (7) 

In this equation Sg is a positive constant, the energy E is measured from the 
Fermi level Ef=0 and q is a real number so that g(E) obeys eq.([2]). Typical 
values of q can be q=l, 3, 5, etc, or q=l/3, 1/5, 1/7 or 5/3, 5/7 etc. 

Let start our conductivity computation by evaluating the number m(Ej), when 
o"jj>o"c- For Ej >0 eq.Q yields: 



where the Ii, I2 and I3 integrals are given by: 

h = J,^^iA-E,fNiE,)dE, 

h = Ji^{A-E,rN{E,)dE, 

h = /-A+i..(A -E^ + E^fN{E,)dE, 

By inserting the DOS of eq.© in Ii, I2 and I3 and by putting x=Ej/A, a 
new dimensionless energy variable, eq.® yields after an extremely laborious 
calculation: 



m(x, g) = mo(x) + g) (9) 

We have written in eq.([9]) the number of conductances attached to the site of 
energy Ej as the sum of two terms: mo(a;) represents the number of conduc- 
tances resulting from the symmetrical part N(Eir) of the density of states and 
ix{x, q) is the number of conductances resulting from the asymmetrical part of 
N(E): 



mo{x) = ^{l + x){l-xf (10) 



Mx,g) = Mo^^^^^/(x,g) (li; 
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Here -Eq is an energy below Ep (see Fig{T]), solution of the equation N(E) = 0; 
r(z) is the gamma function of z and the parameters Mq and f{x, q) are given 
by the following expressions: 



f{x, q) = l- ^ 5^.a;g+3 _ - (1 - 

where we have: Ag = (g^ + 7q + 12) /2; Bg = q^ + 6q + 8; Cg = {q^ + 5g + 6) /2 

We are now ready to solve our percolation-conductivity problem by replacing 
eq.([7]) and eq.Q in eq.©: 



In eq. fll2p we have omitted to write the asymmetrical functions since their 
integrations cancel over positive and negative energies. 

Before tackling the tedious resolution of this equation, may be it would be 
advantageous to examine and compare first the terms of each integrand. If 
we place ourselves in the situation shown in figUl where the temperature of 
the VRH system is such as the active energy layer A < E'q > it becomes then 
apparent that we can neglect in eq. fll2p all terms containing the q) number 
since the surface delimited by these terms is negligible in comparison with the 
one corresponding to mKx) when q> 1/3 (see Figj2]). 



Fig. 1. Density of states representation with an active energy layer A < Eq. 



In such a situation, eq. (fT2|) becomes quasi similar to the one found in [3] for 
a constant DOS: 



[mg(x) + 2mo(a;)/i(x, g).z/q.x^ + fi'^{x, q)] dx 



c = 



(12) 
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Fig. 2. Representation of two integrand terms of eq. ()12p in unit of Mg for q=l/3. 
The inset figure shows how much q) in eq. (fT2]) is small. 



mQ{x)dx 



mo{x)dx 



(13) 



Integrating eq. (fT3|) and taking as in [3j c = 1.7, we obtain: 



(14) 



kT \T 

where = 1.8/{k.N{EF).^^Y^^.The numerical value 1.8 of Tg^^ is similar 
to the one found earlier by Pollak in |3]. Consequently, equations and f|T^ 
show that Mott T~^/^ conductivity is valid with a very good approximation 
for any DOS of eq.©, when q> 1/3. 

In the limit q^O (step like DOS at Ep), we lose the approximation that gave 
eq. f|T3|) since mg(x) and 2mo{x)fi{x,0) become comparable. The re-evaluation 
of eq. ffT^ in this case yields: 



[mKx) + 2mg(x)/i(x, 0).i^o] dx 



[mo{x) + fj,{x, 0).i^o] dx 



-1 



(15) 
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Here again the conductivity follows Mott law: a = exp{ — (Tq/T)"^/^} in which 
the new numerical constant of Tg ^''^ takes the value: 1.76. Such a negligible 
change of Tq^^^ is of course indetectable during any experimental mesurement, 
implying the validity of eg. (1141) for any q. 



4 The unified thermopower 



To compute the Seebeck coefficient S, we make use of the formula: 



here e is the electronic charge and 11 is the Peltier heat. Equation [16] is one of 
the classical Kelvin relations of thermoelectricity. For metals it gives [21j : 



3 e 



d\na{E) 
dE 



(17) 



J E=Ef 



In the variable range hopping regime, 11 may be interpreted [5j as the mean 
energy of any site that belongs to the continuous critical path of conductances. 
We can then write it under the form: 



_ J E.N{E).p{E)dE 
~ JN{E).p{E)dE 



The Peltier heat in this expression appears as an energy E, weighted the by a 
probability factor p(E) that is assumed proportional to the number m(E) of 
conductances attached to a site located at energy E. Inserting equations ([7]) 
and ([9]) in eq. lITS!) we obtain: 



^ / [N{EF).E.f,{E, q) + s,.m,{E).E^+^] dE 

J[N{EF).mo{E) + Sg.E'}.fi{E,q)]dE ^ ' 

In this expression, we have omitted again to write the asymmetrical factors 
of E since their integrations cancel over positive and negative energies. For 
q> 1/3 we can also neglect the second denominator term since the active 
layer of states sohcited for the conduction process must verify eqs (1) and (2), 
in which case we have Sq.E'^ < N{Ef). Integrating and inserting the Peltier 
heat of eq. (fT9ll in the thermopower expression of eq. (fT6i) . we obtain: 
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d\nN{E) 



(20) 



E=Ef 



where F{q) = 60 [r(g + l)/r(g + 7)] (g + l)(g + 4). Replacing the active layer 
energy A by its expression of eq. (fT4l) . we find that the T dependence of the 
variable range hopping thermopower is a class given by: 



5 = F(g)^fc%^T^ 

e 



d\nN{E) 



E=Ef 



(21) 



In the case of a linear asymmetry of density of states (q = 1), the classical 
variable range hopping thermopower formula [4-6] is recovered: 



S = F{l)-k{ToT) 



d\nN{E) 
dE 



(22) 



E=Ef 



with the exact numerical factor found in [5]: F(l) 



In the limit oo (constant DOS near Ef), we have S=0 since F(oo)=0. This 
is in agreement with all thermopower theories indicating that S=0 when the 
density of states is symmetric with respect to Fermi level. 

Eq. (1211) is in fact an unified formulation of the thermopower corresponding to 
Mott conductivity since it includes according to q, the following thermopower 
class/behaviors: 



4-1 Integer metallic class 



It is a class obtained from eq. ipT]) for any odd integer q. S~ T^/^ (q= 1) is 
only one element of this class. 



4-2 Fractional metallic class 



It is a class obtained from eq. (l21l) for any fractional q>l/3, satisfying eqs.([T] 
and [2]). We find in this class for q=5/3, the first VRH thermopower, predicted 
in 1969 by Cutler and Mott 
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4-3 Semiconductor class 



It is a class obtained from eq. (l2T|) for any fractional ) with 0<q<l/3. If q 
tends to zero as 1/p, where p is an infinite odd number (step like DOS at Ep), 
we find , the thermopower due to small polarons found earlier by Triberis and 
Friedman [23] for constant DOS band above Ep. 



4-4 The constant thermopower 



The pivotal value q=l/3, delimits the metallic and semiconductor classes. Its 
thermopower is T-independent and is given by: 



S = F{^)^{kToy 



d\nN{E) 



(23) 



J E=Ef 



and replacing Tq 



Ef 



Putting in this expression d\nN(E)/dE^^^ 

by its expression of eq. (fT4l) . we obtain a T-independent formula of the ther- 
mopower: 



S = Pj- (24) 



where /? is a numerical factor of order 1 and : 



is a characteristic length defined by the nature of the DOS. In amorphous 
semiconductors where a constant thermopower of magnitude ~ lOO/xV-ft'^^can 
be observed [T^, we expect a length of order ^ for L. 



5 Conclusion 



We have shown in this work that the thermopower in the VRH regime may 
have according to the DOS shape of eq.([7]) 3 different behaviours: it can 
have a metallic behaviour if q>l/3, it can have a semiconductor behaviour if 
0<q<l/3, or it can be T-independent if q=l/3. 
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